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Abstract 

FuUerenes and nanotubes consist of a large number of carbon atoms sitting on the sites of a 
regular lattice. For pratical reasons it is often useful to approximate the equations on this lattice 
in terms of the continuous equation. At the moment, the best candidate for such an equation 
is the modified non-linear Schrodinger equation. In this paper, we study the modified non-linear 
Schrodinger equation, which arises as continuous equation in a system describing an excitation on 
a hexagonal lattice. This latter system can e.g. describe electron-phonon interaction on fuUerene 
related structures such as the buckminster fullerene and nanotubes. Solutions of this modified non- 
linear Schrodinger equation, which have solitonic character, can be constructed using an Ansatz 
for the wave function, which has time and azimuthal angle dependence introduced previously in 
the study of spinning boson stars and Q-balls. We study these solutions on a sphere, a disc and 
on a cylinder. 

Our construction suggests that non-spinning as well as spinning solutions, which have a wave 
function with an arbitrary number of nodes, exist. We find that this property is closely related 
to the series of well-known mathematical functions, namely the Legendre functions when studying 
the equation on a sphere, the Bessel functions when studying the equation on a disc and the 
trigonometric functions, respectively, when studying the equation on a cylinder. 

PACS numbers: 73.61.Wp, 11.27.-|-d 
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I. INTRODUCTION 



Fullerenes and fuUerene related structures have gained a lot of interest in the past 20 years. 
Fullerenes were discovered for the first time in 1985 and are carbon-cage molecules with a 
large number n of carbon (C) atoms bonded in a nearly spherically symmetric configuration. 
The configuration is such that three of the valence electrons of the carbon atom form the 
bounds with the neighbouring carbon atoms. This means, that one "free" valence electron 
can "hop" along the the n positions in the fuUerene. 

The Cgo fuUerene, the so-called "buckminster" fuUerene, has 60 carbon atoms arranged 
on a lattice with 20 hexagons and 12 pentagons. It has a diameter of (igo ~ 7 A. In general, 
fullerenes with n carbon atoms consists of lattices with 12 pentagons and (| — 10) hexagons, 
where n > 24 has to be even. 

FuUerides are alkali-doped fullerenes, i.e. fullerenes on which alkali metal atoms such 
as rubidium (Rb), potassium (K), caesium (Cs) or sodium (Na) are put. These alkali 
metals donate one electron each. It was observed that such "doping" leads to a metallic or 
even superconducting behaviour j^. The interesting thing about this is that the transition 
temperature Tc of the fuUerides is quite high for superconductors, e.g. Tc = 33 K for a 
RbCs2C6o 01 • In an attempt to explain this high transition temperature, it was found that 
phonon-electron interactions in the fuUerides can explain this phenomenon quite well ^. 

Nanotubes are fullerene related structures, which were discovered for the first time in 
1991 a. TKe.e molecule. Kave cyUndnc^ ...pe t.e .„bo„ .on,, a. a..a„ged on a 
hexagonal grid. The ends of the nanotubes are typically closed by caps of pentagonal rings. 
WhUe the tubes discovered in 1991 were mult i- layered, single-walled tubes were synthesised 
in 1993 0]. The tubes are distinguished by their chirality, diameter and lengths. They have 
average diameters of 1.2 to 1.4 nm and can be up to 2 mm long. Concerning their chirality, 
the tubes are typically put into three different classes: armchair, zigzag or chiral tubes [7[. 
Depending on their structure, the tubes have different mechanical, thermal, optical and 
electrical properties. The chirality determines e.g. whether the tubes behave metallic or 
semiconducting. These facts also lead to the conclusion that a distortion of the hexagonal 
lattice would influence the energy-band gap. A distortion can be achieved in two ways: either 
by external, mechanical forces such as bending, stretching or twisting or by an interaction of 
the lattice with an internal excitation. Trying to explain the dispersion-free energy transport 
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in biopolymers, Davydov realised [8| that the interaction of an internal excitation with the 
distorted lattice (whose distortion was initially caused by the internal excitation) leads to a 
localised state, a soliton jol. For the case of biopolymers, the internal excitation is an amide 
1 vibration. 

Motivated by the research on fuUerene related structures, the interaction of a 2- 
dimensional hexagonal lattice with an excitation caused by an "excess electron" was studied 
recently The Hamiltonian is of Frohlich-type and describes "electron-phonon" inter- 
actions. It was found that the existence of soliton-like structures depends crucially on the 
value of the "electron-phonon" coupling. 

Interestingly, in the continuum limit, this latter system of equations reduces to a modified 
non-linear Schrodinger equation, where the extra term appearing due to the discreteness of 
the lattice and the "electron-phonon" interaction, leads to the stabilisation of the soliton. 
The appearance of this type of equation has been observed also for a similar system on a 
quadratic lattice Consequently, the modified non-linear Schrodinger equation was 



studied on a 2-dimensional plane in and on a sphere in using an Ansatz previously 
introduced in the study of boson stars and Q-balls Non-spinning as well as spinning 
solutions with and without nodes were constructed. 

In this paper, we study the modified non-linear Schrodinger equation arising as contin- 
uum limit of the equations describing the interaction of an electron-like excitation with a 
he xago nal lattice. In Section II, we review some generalities. Using the Ansatz introduced 
in I15I . Ilfil l. we construct solutions for three different cases, which we describe in Sections 

n 

III-V. In Section III, we discuss the results for the solutions on the sphere obtained in [1^ 
from another point of view, namely we give the spectrum of the equation. In Section IV, 
we discuss the solutions on a 2-dimensional disc with fixed radius and demonstrate how the 
solutions obtained on the plane can be understood in this context. In Section V, we dis- 
cuss the equation on a cylinder with fixed radius and fixed length. We give our conclusions 
in Section VI. 
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II. THE MODIFIED NON-LINEAR SCHRODINGER EQUATION: GENERALI- 
TIES 



The starting point is a coupled systems of discrete equations describing the interaction 
of an excitation (e.g. an "excess electron") with a regular lattice in two dimensions. The 
excitation is given in terms of a complex valued scalar field ip, while the lattice displacements 
are described by two fields, one standing for the displacement in x, the other for that in 



11 



121 and 



y direction. These discrete equations were studied for the quadratic lattice in 
for the hexagonal lattice in Remarkably, it was found for the hexagonal lattice that 
in the stationary limit it is possible to replace the full system of equation, in which the 
excitation field ip and the displacement fields are coupled, by a modified discrete non-linear 
Schrodinger equation. This is not possible for the quadratic lattice. 

In the continuum limit, both systems of equations reduce to a modified non-linear 
Schrodinger equation. The new term appearing in this equation, which results from both 
the discreteness of the lattice as well as from the interaction between the lattice and the 
complex valued scalar field ip ("lattice-excitation interaction"), leads to the stabilisation of 
the solution. The difference between the quadratic and the hexagonal case is just the coef- 
ficient in front of the non-linear term in the equation, which, however is of the same order 
of magnitude for the two cases. 

Since we are dealing with fullerene related structures here, we consider the modified non- 
line. S^^^ equation, which an., a. »..uu. 1... . .Ke he.a.ona, ca.e QQ. 
In dimensionless variables it reads: 



t^ + Aij + Agij + g AIV^I^J = . (1) 

g is the coupling constant of the system and determines the "strength" of the non-linear 
character of the equation. Especially, as can be seen from the derivation of the modified 
non-linear Schrodinger equation from the discrete equations the mass of the carbon 
atom is encoded into g. 

III. EQUATION ON THE SPHERE 

In this section, we will discuss the equation above on a 2-dimensional sphere of radius 
Rs aiming to describe "electronic" excitations on fuUerenes and fuUerides. 



The Laplacian operator in this case reads: 

^spKere = + + ^ J ' (2) 

where 6' e [0 : tt] and <y9 e [0 : 27r]. Note that ip = iplt, 6, (f) is a complex valued scalar field. 

Solutions of (P) considered on the sphere with radius Rg can be characterised by their 
norm: 

27V TT 

= j J \i^\^sm9d9d^ (3) 



as well as by their energy: 

27V TV 



E = J J (\^^\'-2g\^\' + \giW\^p\ry^neded^. (4) 





We will construct normalised solutions, i.e. r]"^ = 1. Note that in a different strategy 
was employed. By rescaling the wave-function ip —>■ ^/'/y^? the coupling constant g can be 
eliminated from the equation. (Of course, the normalisation of the wave-function is lost in 
this approach, but since rj was plotted, the corresponding value of g could be easily deter- 
mined.) In this paper, we put the emphasis on the evolution of the spectrum (corresponding 
to normalised eigenfunctions) as a function of the coupling constant g. Note that g deter- 
mines the "strength" of the non-linear part of the equation and that for g = 0, the equation 
becomes the ordinary Schrodinger equation on the sphere, which, of course, is linear. 

A. Discussion of the solutions 

To construct explicit solutions, we employ the following axially symmetric Ansatz jl^ Q : 

^(t,^,</?) = e-*"*+*'"'^$(^) . (5) 

The solutions with m = are non-spinning solutions, while the spinning solutions will have 
m 7^ 0. 

Inserting this Ansatz into (^, we find the following equation [l^ : 

1 + 1 + 5($2 1 + ^($2 ^ ^ 

where the prime denotes the derivative with respect to 6. 
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This is a priori a non-linear equation with external parameters g and u. Here, we will 
consider it from a slightly different point of view. Let us rewrite © as follows: 

^ f*" + cot »*' + *"^ + £^^i<l^- ^-1^*1 =M (7) 



1 + 5«l>2 l + g^'^sin 6 
where 

Q = ujRl (8) 

is the effective spectral parameter. The equation then appears as an eigenvalue Schodinger 
equation supplemented by a non-linear term coupled with the parameter g. The solutions 
in the linear case, i.e. for (7 = are well-known and it is very likely that - by continuity- 
their deformation due to the non-linear term can be constructed at least for small values 
of g. However, this does not exclude the occurence of supplementary branches of solutions, 
which may appear typically due to the non-linear structure of the equation. 

In order to construct regular solutions, (jT)) has to be solved subject to the following 
boundary conditions: 

$1e=o = $1e=. = for m = (9) 

and 

$(0 = 0) = $(0 = tt) = for m ^ . (10) 

To solve the equation ((7j) numerically subject to the boundary conditions 0, respectively 
(I10|l for generic values of g and i?^, we use a numerical routine based on the Newton- Raphson 
method p/O]. Our strategy is to use $(6' = 0), respectively ^'\e=Q as a "shooting" parameter 
for the case m = and m 7^ 0. 

1. The case g = 

In the case g = 0, ^ reduces to the well-known Legendre equation 



<l>" + cot^$'-m2^^$ + ^]$ = . (11) 
sin 6 



The solutions of this equation are the associated Legendre functions: 

$(^) = Pl^{cose) with Q = l{l + 1) , (12) 

where 

I EN, m = -/,-/ + 1,...,/- 1,/ . (13) 



The spatial part of ^ are then the spherical harmonics Y^{9,ip) = e*"*'^P^(cos^). Notice, 
however, that in the equation studied here there is no distinction between m and —m. 
From now on, we will assume m to be positive, but it is understood that the solution with 
corresponding —m can be trivially obtained. 

The spherical harmonics have a well-defined parity, i.e. they are either symmetric or 
anti-symmetric under the refiection ^{6) ^{tc — 6). In addition, the number of zeros in 
the interval ]0 : 7r[ is equal to / — \m\. 

These results are, of course, very well- known, but we mention them for completeness 
and because the spectrum of the full equation will approach the spectrum of the Legendre 
equation in the g ^ limit. 

2. Deformed Legendre functions 

All the solutions discussed in the previous section are continuously deformed when the 
parameter g is chosen to be larger than zero. For every available spherical harmonic, we 
observe that a branch of solutions in g exists, on which the solutions have the same symmetry 
than in the g = limit. Our numerical results further indicate that the number of zeros k 
of the solution is preserved all along the branch and is equal to k = I — m. 

Now, we will discuss in detail the numerical solutions corresponding to the deformations 
of the spherical harmonics Yq, and Y^. 

• The non-spinning (m = 0) and nodeless {I = 0, k = 0) solution corresponding to the 
case Yq is in fact the constant solution of the equation, namely 



The eigenvalue f2 of the corresponding normalised solution has a linear dependence of 
the parameter g: 



This dependence is shown in Fig^for Rs = 2. 

• The non-spinning (m = 0), one-node [l = 1, k = 1) solution is anti-symmetric under 
the refiection 6 ^ n — 6. The function Yq oc cos 6' is progressively deformed by the 
non-linear term for g > 0. For g ^ 10 the profile of the function $(6') possesses 




(14) 
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a plateau $ = surrounding the region of 9 = tt/2. Away from this plateau, the 
solutions resemble two disconnected solitons located at the two poles of the sphere, 
i.e. the maxima of $^(^) are situated around 6 = n and 6 = 0. Stated in other 
words, the probability density vanishes in a region surrounding the equator, while it 
has local maxima at the two poles. Let us stress, that the physical parameters remain 
monotonic functions of the "shooting" parameter = 0). 

The dependence of Q on g is shown in FigQfor Rg = 2. For g = 0, we have Q = 2, as 
expected. It is worth noticing that the two eigenvalues corresponding to the quantum 
numbers m = I = and m = 0, / = 1 cross at ~ 1.5. This would a priori mean, 
that the "excited" solution becomes the "ground state" of the equation. However, as 
we will see in the next section, the appearance of new asymmetric solutions will turn 
out to have even lower values of Q. 

• The spinning (m = 1), no-node (/ = 1, A; = 0) solutions arrive as deformations of the 
spherical harmonic Y"/ oc sin^e*'^. This solution is thus symmetric under the reflection 
6 ^ 71 — 9. In contrast to the non-spinning solutions, we notice here that when 
the coupling constant g is of order 2.5, the structure of the solutions becomes very 
sensitive to the value of the "shooting" parameter $'|6i=o- In other words, there exists 
a critical value of ^'\e=o = $(„.)|6i=o (depending on R^) such that no solutions exist for 
$'|e=o > $(cr)|0=o- Our numerical results, however, suggest that solutions of arbitrary 
large values of g can be constructed when considering the limit ^'\e=o ^[cr)\s=o- The 
corresponding eigenvalue of Q decreases monotonically with g from Q = 2 (for the case 
g = 0) as can be seen in Fig^for the case Rg = 2. Of course, this curve also represents 
the solution with m = —1, so that the eigenvalue fl is three- fold degenerated in the 
limit g = 0. 

Finally, let us state that for these solutions, the function $ has a maximum at the 
equator and vanishes at the poles. This is an interesting result with view to 18]. In 
this latter paper it was found that if a Ceo or a C70 molecule is doped with one or two 
excess electrons, the additional charges accumulate nearly along an equatorial line of 
the molecule. 

Of course, deformations of the Legendre functions for higher m and / can be constructed 
in a systematic way. 
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3. Asymmetric solutions 

All solutions of the non-linear equation (|7j) discussed previously appear as continuous 
deformations of the solutions available in the linear limit g = 0. However, in this section, 
we will exhibit new branches of solutions, which are fundamentally non-linear phenomena 
and as these have no linear counterparts. Numerically, we find that these new branches 
of solutions exist for sufficiently high values of the coupling constant g parametrising the 
non-linear effect and that some of these branches bifurcate into the branches of the deformed 
Legendre functions at critical values of g = gi^'''\Rs), which depend on I, m (and thus on 
the number of nodes k = I — m) as well as on the radius Rg. 

These solutions are charaterised by the fact that they are asymmetric with respect to the 
reflection 9 ^ n — 9. 

• Asymmetric, non-spinning (m = 0) and nodeless (/ = 0, A; = 0) solutions exist for 
g > gcT^\Rs = 2) ^ 0.838. When the parameter g is chosen larger than this critical 
value, our numerical results indicate that this solution has a maximum at ^ = and 
a minimum at = vr and decreases monotonically between these two extrema, which 
are strictly positive. Thus, the solution is neither symmetric nor anti-symmetric under 
the reflection 9 ^ it — 9. Instead, this reflection applied to the solution leads to an 
equivalent solution with a maximum (respectively minimum) at 6' = vr {9 = 0). The 
solutions thus represent a physical state with probability density concentrated in one 
of the hemispheres. This is a new feature that cannot arise in the deformation of the 
Legendre functions. 

In the limit g —>■ gcr^^ the values of the maximum and the minimum approach each 
other and the solution becomes constant. The branch of these asymmetric solutions 
bifurcates into the branch of constant solutions (deformed Yq) described in the previ- 
ous section. The corresponding value of u of these asymmetric solutions is lower than 
the one of the symmetric branch. This is shown in FigG] (dotted line) for Rg = 2. 

For generic values of Rs, the bifurcation of the asymmetric branch occurs at 




71 



(16) 




with 





(17) 



R's - 1/4 ■ 
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• Asymmetric, non-spinning (m = 0) and one-node (/ = 1, = 1) solutions exist for 
g > gcr''^\Rs = 2) ^ 25. In this region of parameter values, we find that in fact two 
branches of solutions exist (see Fig^. For a fixed value of g, the solution with the 
lower value of Q has a bigger value of $(0) and its minimum lies at smaller values of 
6. This is demonstrated for g = 37, i?^ = 2 in Fig|21 The values of u corresponding to 
these asymmetric solutions are higher than the values of the corresponding (i.e. with 
for the same values of g) antisymmetric solution on the branch of deformed Legendre 
functions Yq. In contrast to the other described cases, we have not succeeded to 
construct another (a third) branch which would eventually bifurcate from the branch 
of deformed Legendre functions. 

• Asymmetric, spinning (m = 1) and nodeless (/ = 1, = 0) solutions exist for 
g > gcT^\Rs = 2) ^ 3.79. Solving the equations for values of g larger than this 
critical value, we were able to construct a branch of asymmetric solutions (and the 
corresponding mirror symmetric one). The maximum of these solutions is located be- 
tween ^ = and 9 = 71/2 such that the probability density is located in one of the 
hemispheres. As before, the frequency of the asymmetric solution is lower than the 
frequency of the corresponding (i.e. same g) solution on the deformed branch. 

The construction of asymmetric branches with higher values of the spin m and/or higher 
node numbers k becomes very tricky. This is mainly due to the occurence of many solutions 
existing for the same values of the parameters and also due to the fact that the profiles of 
the symmetric and asymmetric (or anti-symmetric and asymmetric) solutions can be very 
close to each other on a large part of the [0 : vr] interval. The construction of asymmetric 
branches can thus not be achieved systematically. 



IV. EQUATION ON THE DISC 

In the equation was discussed on a two-dimensional plane and non-spinning as 
well as spinning solutions with nodes were constructed. Here, we will treat this problem 
differently. With view to the discussion of the equation on a sphere, we will discuss here the 
equation on a two-dimensional disc with radius Rd parametrized by the coordinates p and 
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lied in 



yj. The case studied in jl3| then arises as hmit Rd oo. The Laplacian is given by: 

. _ Id 1 

^disc — TT^ H — t: — I — ^'TT^ (18) 
Op^ p op p'^ Oif^ 



with G [0 : 2tx\. The norm reads: 



27r Rd 

= J J m'pdpdf (19) 





and the energy is given by: 

E = j j + \9{^\i^?f^ Pdpd^ . (20) 



A. Discussion of the solutions 

As for the case on the sphere, we adopt the following Ansatz for the complex valued 
function %lj{t, p, ip): 

^{t,p,^) = e-'^'+'"^^<P{p) . (21) 
Inserting this Ansatz gives the following equation: 

where the prime now denotes the derivative with respect to p. 

We can rewrite, analogue to before, the equation in the following way: 

This equation looks like the one studied in however, we impose different boundary 
conditions here: 

9p0|,=o = , 0(i?rf) = for m = (24) 

and 

0(p = 0) = , (f){Rd) = for m ^ . (25) 

Again, we used the numerical routine described in [l7| to solve the equation subject to 
the boundary conditions (f^ . respectively (|^. 
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1. The g = limit 



In the limit g = (PHj) reduces to the hnear Schrodinger equation in polar coordinates 

0" + ^-^ + -0 = O. (26) 

Note that for static solutions, i.e. u = 0, the equation reduces to the Euler equation with 
solutions of the form ip{t^p,ip) = ipip) oc e^'^'^p"*. For u ^ 0, the equation becomes the 
Bessel equation if we define x = \/\uj\p and F{x) = (j){x): 



cPF IdF f 



dx^ ^ X dx^ \ x'^j ^ ^ 

The solutions of this equation, which are regular at the origin are of course the well-known 
Bessel functions Jm{x) for uj > and the modified Bessel functions Im{x) for uj < 0. 
Remembering that the Bessel functions Jm{x) are oscillating and admit an infinite number 
of nodes on the positive real axis, it is easy to construct solutions of with k nodes. This 
is possible assuming that 

0(p) = JmiV^p) with = (28) 

and choosing the value of u in such the way that the value p = Rd coincides with the /c-th 
zero (we call it x = Xm,k ) of the Bessel function J^, i.e. y/u = Xm,k/Rd- 

This indicates that in the limit g = ^ solutions with angular momentum m and with k 
nodes exist and the spectrum of the equation is positive. The corresponding functions are 
convex. 



2. The g ^ case 

If the coupling constant g is positive but small, we expect that the pattern of solutions 
will be similar to the one available in the limit g = discussed above. Our numerical results 
indeed confirm this expectation. For small values of g, the solutions are still represented by 
convex functions and the values of u are positive and decrease slowly with g as shown in Fig. 
EJfor the case m = 0, k = 0. Let is call this range of (7- values the "weak coupling-regime". 

When g is large enough, the value of u becomes negative (see FigO)) and the solutions 
have an inflexion point at some intermediate value of p (see Fig0)). As a consequence they 
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become much more localized in the region around the origin. In this "strong coupling- 
regime" the solutions can thus be described as being soliton-like. We also note that the 
value of oj decreases very strongly with g. 

The difference between the strong and weak coupling is shown in FigEJfor (7 = 1, respec- 
tively g = 8 and Rd = 5. The value u is shown in dependence on g for the non- spinning 
(m = 0) and nodeless {k = 0) solution for different values of the disc-radius Rd in FigOl 
In Fig. Elwe also superposed the values of uj corresponding to the case of the full plane (i.e. 



Rd = oo). As has been shown 



in |lO 



13l |. in this case, solutions exists only for large enough 
values of g: g > gcr ~ 2.94. It is therefore natural to try to understand the pattern available 
on the plane as the Rd oo limit of the pattern of solutions for the disc. In this respect, it is 
interesting to notice that in the "strong coupling-regime" the values of u corresponding to a 
finite disc get very close (especially for the case corresponding to Rd = 10) to the values of uj 
associated with the case of the infinite disc. This is related to the fact that with increasing 
g, the "soliton" becomes more and more peaked around the origin and doesn't "notice" 
whether the disc is finite or infinite. This observation gives a natural explanation for the 
fact that solutions exist only for g > gcr in the case of the plane. Indeed in the case of the 
plane, it was noticed that the linearized equation (valid asymptotically) also is of Bessel 
type. However, the nodeless, localized solutions are exponentially decreasing and can only 
be associated with the modified Bessel function, which correspond to negative values of u. 
That's why, on the plane, nodeless soliton-like solutions exist only for g > 2.94 and have no 
"weak coupling" -counterparts. For large values of g the numerical results strongly suggest 
the following linear relation between g and the frequency corresponding to the fundamental 
mode : u ^ —0.66g + 2.19. We also studied the solutions with k > nodes and m > and 
found that the pattern is similar to the case k = and m = discussed above in detail. 

We do not find it useful to discuss them here, but we believe that our results provide 
a robust evidence that regular solutions of p3|l with an arbitrary number of nodes and 
arbitrary (integer) angular momentum exist on the disc/plane. 

V. EQUATION ON THE CYLINDER 

We also considered the equation on a cylinder of radius Rc and length 2L. This is a natural 
extension of the previous results in order to describe soliton-like structures on nanotubes 
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with radius Rc and length 2L. Using the natural coordinates z, (p, the Laplacian reads: 

1 92 Q2 

^cylinder = ^^9^+9^ (29) 

with ip E [0 : 27r] and z G [— L, L]. The norm reads: 

2tv L 

= J j l^l^dzdif (30) 

-L 

and the energy is given by: 

271 L 

E = j j {\Vi^\'-2gW + \g{V\i^\'Y^dzd^. (31) 

-L 

A. Discussion of the solutions 

We use the following Ansatz: 

V'(t, 2, ^) = e-*^*+*™'^F(2) , m = 0,1,2,... (32) 
The equation (jT)) then becomes: 



^" + (^r^+^49f^ + --^j^=0. (33) 

where the prime denotes the derivative with respect to z. 

Again, we can rewrite this equation in terms of a "spectral equation" : 



coF . (34) 



This equation has to be solved subject to the following boundary conditions: 

F{-L) = F{L) = . (35) 
The equation is then solved subject to the boundary conditions using the numerical 



routine described in 



17|. 



1. The case g = 



The case g = corresponds to an harmonic equation: 

2 

m 



F" = -(.--)F. (36) 
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The solutions of this equation subject to the boundary conditions are: 
F(z)=sin((z-L)v/^-m7i?2) with ^ = + (^l)' ' ^=1'2,3,... (37) 

B. The case 5 / 

It is possible to solve by quadrature, which gives: 

and C is an integration constant. The integration can be performed, however, it leads to 
an involved and untractable expression of z in terms of F, C, Q. Assuming the cylinder to 
be infinite {L = 00) and F{oo) = F'(oo) = leads to C = 0. If we now require in addition 
that the solution is symmetric under the reflexion z —z (i.e. F'{0) = 0) we obtain the 
relation Q = —2gF{0y. This provides a useful check for our numerical results. 

Performing a numerical analysis of the values of as a function of g (with L, R^, k, m 
fixed), we obtain a pattern very similar to that available for the case on a disc (see FigE)). The 
only noticable difference resides in the fact that, in the limit L = 00, a normalizable solution 
can be constructed for all values of the parameter g. This case is, of course, interesting with 
iew to nanotubes since the lengths 2L of a tube is much larger than its diameter 2Rc. In the 
case of the "ground state" solution, which is non-spinning (m = 0) and nodeless {k = 0), 
the following relations hold : 



00 



-O.lg^ for ^ < 1 , -0.555c/ + 1.28 for ^ > 1 . (39) 



Unlike in the case of the sphere, we did not succeed in constructing asymmetric solutions in 
the case of the cylinder (i.e. solutions, which are neither even nor odd under the reflexion 
z —>■ —z) and we believe that there are no such solutions. 



(40) 



VI. CONCLUSIONS 



In this paper we have presented an extended analysis of a modifled non-linear Schrodinger 
equation in 1+2 dimensions. This equation was constructed as an effective continuous 
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equation describing an interaction of an excitation with a 2-dimensional hexagonal lattice. In 
this paper, we think of the excitation as an "excess electron" described by a wave function ip 
and the interaction to be "electron-phonon" -like. As is well-known from amide I-excitations 
in biopolymeric lattices such an interaction leads to the creation of a localised state, a 
"soliton" . 

The continuous modified non-linear Schrodinger equation was obtained with the idea 



10, ll^, but it can be modified in order 



to approximate the physical system on the plane 
to be considered on different geometrical manifolds, which are justified physically. These 
are, namely, a cylinder to mimick excitations on a carbon nanotube and a sphere, respec- 
tively, in order to describe excitations (e.g. extra electrons) on a fuUerene nanomolecule, 
which is particularly interesting with view to the possibility of describing the high transition 
temperature of superconducting fullerides. 

Completed with the different types of boundary conditions related to the geometry and 
using an Ansatz previously introduced m the study of boson stars Q and spinnning Q-balls 
|lfil |. the equation ((H) possesses a lot symmetries, some are continuous (like rotations), some 
are discrete (like reflexions). 

In all cases considered the spectrum of solutions is rich, many "normal-mode" -types of 
solutions exist, mainly characterized by two integers: the internal angular momentum (or 
spin) and the number of nodes of the wave function. This is very reminiscent to the principal 
quantum numbers of more familiar quantum mechanical systems. 

It might turn out to be difficult to classify the solutions of a non-linear, partial differential 
equations. The various types of solutions, however, can be traced back from the correspond- 
ing solutions available in the weak coupling limit, in which the equation becomes linear and 
its spectrum is known analytically. In the case of the sphere, disc and cylinder, respectively, 
we found deformed Legendre functions, deformed Bessel functions and deformed harmonic 
functions. However, in the case of the sphere, we demonstrated the existence of additional 
solutions, which are specifically related to the non-linear character of the equation. The 
main feature distinguishing these supplementary solutions is their asymmetry with respect 
to the refiexion at the equator of the sphere. In other words, we exhibited the spontaneous 
symmetry breaking (SSB) of one of the natural symmetries of the initial problem. 

This SSB occurs when the non-linear coupling becomes strong enough. Similarly it can 
occur for very different types of non-linear equations. As an example, we point out the 
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bifurcation of bisphaleron solutions from the sphaleron solutions available in the case of the 
classical equations of the electroweak field theory 20 1. In this case, the parity operator 
is spontaneously broken by the so called bisphaleron solutions. 
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deformed Legendre 




1 2 3 4 5 6 



g 



FIG. 1: The eigenvalue = coRg is shown as a function of g for {m,l) = (0,0), (0, 1) and (1,0). 
We have chosen Rg = 2. The sohd Unes denote the deformed Legendre functions, while the dotted 
lines are the asymmetric solutions (see text for more details). 
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FIG. 2: The profiles of tlie function ^(6) and the derivative ^ are shown for {m,l) = (0,1) 
(non-spinning, nodeless solutions), g = 37, sphere-radius Rg = 2 and two different values of the 
spectral parameter = uRg. 
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FIG. 3: The eigenvalue u is shown as a function of g for the non-spinning (m = 0), nodeless 
{k = 0) solutions on the disc. We have chosen = 5, 10 and oo. Note that = oo corresponds 
to the case of the plane. 
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FIG. 4: The profiles of tfie function 0(p) and the derivative ^ are shown for non-spinning 
nodeless {k = 0) solutions on the disc with radius Rd = 5 for two different values of g. 
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